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Abstract 

We show how a lattice approach can be used to derive Thermodynamic Bethe Ansatz 
(TBA) equations describing all excitations for boundary flows. The method is illus- 
trated for a prototypical flow of the tricritical Ising model by considering the continuum 
scaling limit of the A4 lattice model with integrable boundaries. Fixing the bulk weights 
to their critical values, the integrable boundary weights admit two boundary fields £ 
and n which play the role of the perturbing boundary fields cpi t 3 and (pi 2 inducing 
the renormalization group flow between boundary fixed points. The excitations are 
completely classified in terms of (m, n) systems and quantum numbers but the string 
content changes by certain mechanisms along the flow. For our prototypical example, 
we identify these mechanisms and the induced map between the relevant unitized Vira- 
soro characters. We also solve the boundary TBA equations numerically to determine 
the flows for the leading excitations. 



1 Introduction 

A problem of much current interest in Quantum Field Theory (QFT) is the renormalization 
group (RG) flow between different boundary fixed points of a Conformal Field Theory (CFT) 
that remains conformal in the bulk. Many boundary flows have been studied including |T], [J 
Lee- Yang, minimal models and c = I CFTs. In addition, approximate numerical scaling 
energies can be explored by use of the Truncated Conformal Space Approach (TCSA) ||. 

In this letter, for boundary flows, we show how one can derive exact TBA equations [[| 
for all excitations from a lattice approach 0, 0, §[• This has not been possible in a field 
theory approach. We illustrate our approach in the context of a prototypical boundary 
flow of the Tricritical Ising Model (TIM) || [T(| with central charge c = 7/10. This is 
an interesting but relatively simple CFT with many applications in statistical mechanics. 
Following Cardy [11], the conformal boundary conditions -B( r , s ) of the TIM are in one-to-one 
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correspondence with the six primary fields <£(»•,«) = V(4-r,5-s) labelled by the Kac labels (r, s) 
with r = 1,2,3 and s = 1,2,3,4. The Kac table of conformal weights A r s is shown in 
Figure |I[ Of the many possible flows on a cylinder between boundary conditions -B( r ,s)|(r',s') 
(that is -B( r ,s) on the left and Bi r i iS /\ on the right) we only consider here flows involving a 
non-trivial boundary Bt ra \ on the left with a fixed trivial boundary -6(1,1) on the right. For 
these conformal boundary conditions on the cylinder we use the special notation 



B 



(r,s) 



B 



In Figure [TJ the possible integrable flows [TIJ between these conformal boundary conditions 
are shown schematically. 
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Figure 1: Kac table of conformal weights A r s = ((5r — 4s) 2 — l)/80 of the tricritical Ising 
model. We restrict to the lower half of the table with s = 1,2 using the usual identification 
(r,s) = (4 — r, 5 — s). The integrable flows between the six conformal boundary conditions Br r)S ) 
are shown schematically. The 95(1,3) and ^(1,2) flows are indicated by solid and dotted arrows 
respectively. We omit non-Cardy type boundary conditions. 

In this letter we will focus on the prototypical flow obtained by perturbing the $(1,2) 
boundary TIM by the relevant boundary operator 92(1,3) 

UV = B m ^ %d = IR, XiM - X2,i(q) (1.1) 

This flow induces a map between Virasoro characters Xr,s{<l) of the theory where q is the 
modular parameter. The physical direction of the flow from the ultraviolet (UV) to the 
infrared (IR) is given by the relevant perturbations and is consistent with the g-theorem 
which asserts that the boundary entropy is reduced under the flow gm < guy. The boundary 
entropies are 



/2\ 1/4 sin^ sin^ r r 

9(r,s) = r ===== , ^(1,1) : </( 2 ,i) : ^1,2)^(2,2) = 1 : V2 : g : V2g (1.2) 

V o j u sin j sin t 



4 u '" 5 



with g = (1 + V / 5)/2 the golden mean. 

A more complete account of the full set of six integrable flows of the TIM, including 
the second flow emanating from the £>(i,2) boundary condition, will be given in a subsequent 
paper |K|. The approach outlined here, however, is quite general and should apply, for 



example, to all integrable boundary flows of minimal models [14 
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2 Scaling of Critical A4 Model with Boundary Fields 



It is well known that the TIM is obtained as the continuum scaling limit of the generalized 
hard square model of Baxter [pj| on the Regime III/IV critical line. This integrable lattice 
model is the RSOS lattice model of Andrews-Baxter- Forrester ||16|| . Specifically, in the 
presence of integrable boundaries, the scaling energies of the TIM are obtained from the 
scaling limit of the eigenvalues of the commuting double-row transfer matrices D{u) [T7J of 



the A4 lattice model with N faces in a row. It has been shown by Behrend and Pearce 
that at criticality there exist integrable boundary conditions associated with each conformal 
boundary condition (r, s) carrying 0,1 or 2 arbitrary boundary fields £, 77. In the present case 
the boundaries (1, 1) and (3, 1) admit no fields, (2, 1) and (1, 2) admit a single field £ and 
(2, 2) admits two boundary fields £ and r\. These boundary fields induce the ^(1,3) an d ¥>(i,2) 
RG flows respectively so we identify the thermal and magnetic boundary fields accordingly 

£~<£>(1,3), »7~¥>(i,2)- (2-1) 
To consider boundary flows we must fix the bulk at its critical point and vary the 
boundary fields £, 77. But paradoxically, it has been shown |J that the fields £, 77 are irrelevant 
in the sense that they do not change the scaling energies if £ and 77 are real and restricted 



to appropriate intervals. Explicitly, the cylinder partition functions obtained || [L8| from 
(r, s) integrable boundaries with £, 77 real are independent of £, 77 and given by single Virasoro 
characters 

B(r, s ) = £(r,s)|(l,l) : Z(r, S )|(l,l)(<?) = Xr,s(<?)- (2.2) 

The reason for this is that in the lattice model the fields £ and 77 control the location of zeros 
on the real axis in the complex plane of the spectral parameter u whereas only the zeros in 
the scaling regime, a distance i log iV out from the real axis, contribute in the scaling limit 
N — > oo. The solution is to scale the imaginary part of £, 77 as log N. This is allowed because 
£, 77 are arbitrary complex fields. 

For our prototypical flow (|1.1|), following || we consider a cylindrical lattice with 
the (1, 1) boundary on the right (with no boundary field) and the (2, 1) boundary on the left 
with boundary field £l- Explicitly, we now scale £^ as 

Re(£zJ = - = - , Im(£i) = (2.3) 

where £ is real. In terms of boundary weights, the (1, 2) conformal boundary is reproduced 
in the limit Im(£^) — > ±00 whereas the (2, 1) conformal boundary is reproduced in the limit 
Im(&)-0. 

The scaling behaviour of the A 4 lattice model at the boundary fixed points is described 
by the TBA equations of O'Brien, Pearce and Warnaar |J (OPW). In this letter, we gen- 
eralize their analysis to include the scaling boundary field ( |2.3| ). This parameter appears in 
the TBA equations only through the boundary contribution g(u) in (6.3) of OPW. For the 
flow fll.lf) , the scaling limit of this function in the first analyticity strip — f < Re(tt) < ^ 
now becomes 

^(x,0 = coth^ (2.4) 
with £ — > —00 and £ — > +00 corresponding to £>(i,2) and #(2,1) respectively. 
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3 Classification and Flow of Excited States B(\ 2) >— ► ^(2,1) 

As in OPW, the scaling energies are classified by the patterns of zeros of the eigenvalues of 
the double row transfer matrix D(u). In the complex u plane, they are organized in two 
analyticity strips — | < Re(u) < ^ an d 2A < Re(ii) < 4A. For N sufficiently large, there 
can be only simple zeros in the middle of each strip (1-strings) or pairs of zeros located 
on the two borders of each strip, having the same imaginary part (2-strings). Because of 
crossing symmetry, all the zeros appear in complex conjugate pairs so the upper and lower 
half planes are related by reflection in the real axis and this remains true for general values 
oflm(a). 

The patterns of zeros are classified || by string content and ordering. Let > be 
the number of 1-strings in the upper half plane in strips 1 and 2 and similarly rii , 2 > the 
number of 2-strings. Then for the boundary fixed points Sri 2) an d $(2,1), the string contents 
are classified by the (m, n) systems 

m/ n ( aa \ N+m 2 + a m x -a 

UV = £5(1,2) (mi odd, m 2 even) : n x = - mi, n 2 = m 2 (3.1) 

IR = £>( 2 ,i) (mi odd, m 2 odd) : n t = N J ^ m ' 2 _ m ^ n2 — !^L-t 1 _ m2 (3.2) 

where N is odd. In the UV there are two distinct (m, n) systems labelled by a = ±1. In the 
scaling limit iV — > 00, we have n\ ~ N/2 with mi,m 2 ,n 2 = 0(1). For given string content, 
the relative ordering of 1- and 2- strings is given by the quantum numbers lj® where, for the 
1-string yj£\ 1^ is the number of two-strings whose location zf 1 is larger (further from the 
real axis) than yj£\ Labelling the strings so that y^ +l > y%\ the quantum numbers must 
satisfy 

n % > I? > I? > ... > J» > 0. (3.3) 

In B(i 7 2) with <T = 1 there is a frozen 1-string, that is, the 1-string furthest from the real axis 
in strip 1 is restricted to have no two-strings above it so that Iml = 0. With these quantum 
numbers, the states are uniquely labelled as 

UV = £(1,2) : / = (A (1) I? ... I A (2) If ... I%1 ). (3.4) 

IR = Bfti) : / = (I? I? - IS 1 1? I? - til ) (3-5) 

In the scaling limit, n\ is infinite so in this limit the corresponding strip 1 quantum numbers 
are unbounded. 

Given the pattern of zeros, that is, the string content and quantum numbers at the 
boundary fixed points, the conformal dimensions or scaling energies of these states are 
uniquely determined M by 

UV = B llA : E = A {l!2) + n uv =± + \rnCrn-^^a + Y: £4° (3-6) 

i=l,2 k=l,..,rrii 

IR = B (2A) : E = A i2A)+ n IR =^- 1 -+ 1 -mCm+J2 £ 4° (3-7) 

1=1,2 fe=l,..,m; 
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Figure 2: The three mechanisms A, B, C respectively that change string content during the flow 
^(2,1) l— ^ ^(1,2)- Note that this is the reverse of the physical flow. These mechanisms are illustrated 
for the states: A, (0|0) ■-> (0)+; B, (1|0) ■-> (0)_; C, (000|1) ■-> (000 |00)_. 



where m = (mi, 7712) and the Cartan matrix is 




From 0, we know that the above counting of states reproduces exactly the expected UV and 
IR characters at the boundary fixed points but there must be some mechanisms to change 
these classifications during the flow. 



3.1 Three mechanisms for changing string content 

For convenience, in this section only, we consider the flow IR = £>(2,i) 6(1,2) — UV which 
is the reverse of the physical flow ( |1 . 1| ) . By simple counting arguments, the eigenvalues of 
the double row transfer matrix have 2N + 4 zeros in the periodicity strip for 6(2,1) but only 
2N + 2 zeros for 6(1,2)- Thus precisely two zeros, one in the upper half plane and one in 
the lower half plane, must "flow off to infinity" during the flow. Indeed, from empirical 
observations based on direct numerical diagonalization of small double row transfer matrices 
we find just three mechanisms for changing the string content during the flow. These all 
involve the farthest zeros (1 or 2 strings) from the real axis in the upper half plane and are 
well supported by numerical analysis of the TBA equations in Section 5. 

The mechanisms for the flow IR = 6(2,1) l— * 6(1,2) = UV, as shown in Figure 0, are: 

A. The top 1-string in strip 2 flows to +00, decoupling from the system while mi, n\ and 
n 2 remain unchanged. This mechanism applies in the IR when Iml = Iml = and 
produces frozen states in the UV with a = 1 and Iml = 0. 

B. The top 2-string in strip 1 and the top 1-string in strip 2 flow to +00 and a 2-string 

comes in from +00 in strip 2 becoming the top 2-string. Consequently, each 7j 

(2) 

decreases by 1 and each 7^ increases by 1. This mechanism applies in the IR when 

Iml > an d Iml = and produces states in the UV with a = — 1 and either 7712 = 
or ig > 0. 

C. The top 2-string in strip 2 flows to +00 and a 1-string in strip 2 comes in from +00. 

Consequently, each 7^ decreases by 1. This mechanism applies in the IR when 7m 2 > 

(2) 

and produces states in the UV with a = — 1 and 7m 2 = 0. 
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Observe that these mappings are in fact one-to-one, that the 1-strings in strip 1 are never 
involved in these mechanisms and that the change in parity of m2 is consistent with ( |3.1| ) 
and (IS - ! 



A,B: 



in 



uv 



u 2 — m 2 — j-, w . no 2 — nv 2 

The explicit mapping of the first 18 states is shown in Table 1. The consistency and com- 
pleteness of these three mechanisms is demonstrated by the explicit mapping between the 
finitized characters associated with the UV and IR fixed points. The situation here is very 
similar to the massless thermal <£>i 3 flow M from the TIM to the critical Ising model and is 



IR 



- 1, 



C : m 



uv 



m 



IB 



+ 1 



(3.9) 



best described in terms of finitized Virasoro characters [19| 



n 


Mapping of states and mechanism 


n 


Mapping of states and mechanism 





(0) + - 


- (0|0) 


A 


5 


(4)- - 


- (5|0) 


B 


1 


(0)- - 


-> (1|0) 


B 


5 


(110)+ h 


-> (1 1 0|0) 


A 


2 


(1)- - 


- (2|0) 


B 


5 


(2 0)+ h 


-> (200|0) 


A 


3 


(2)- - 


- (3|0) 


B 


6 


(110|00)_ h 


-> (1 1 0| 1) 


C 


3 


(0 0)+ 1- 


-> (000|0) 


A 


6 


(200|00)_ h 


-> (200|1) 


C 


4 


(000|00)_ 1- 


-> (000|1) 


C 


6 


(5)_ H 


- (6|0) 


B 


4 


(3)_ i~ 


- (4|0) 


B 


6 


(0 0)_ h 


- (1H|0) 


B 


4 


(10 0)+ H 


-> (100|0) 


A 


6 


(210)+ 1- 


-> (2 1 0|0) 


A 


5 


(100|00)_ h 


-> (100|1) 


C 


6 


(3 0)+ h 


-> (3 0|0) 


A 



Table 1: Explicit mapping of states in the order of increasing energy for the first 18 energy 
levels, where n = n uv = n IR is the excitation level as in ([O]) and ( |3.7[) . For some higher levels 
n uv + n IR . 



3.2 RG mapping between finitized characters 



Finitized characters are the generating functions for the scaling energies of the A4 model for 
finite M x N lattices: 

M/2 



Ar.s 



w = ^ c/24 E^ = ^ c/24Tr 



D (N) { 



(3.10) 



where the modular parameter is q = exp (— tx sin(5w) M/iV) and D^ N \u) is the largest eigen- 
value of the double-row transfer matrix D^(tt) with (r,s) boundary conditions. The fini- 
tized characters involve Gaussian polynomials which satisfy a recursion used repeatedly 
below 

n Ii Im-l 



m+n 

m \ 



h =0/2=0 lm=0 



+In 



rn- 



n-l 
m—1 



+ q m [ 



m rn— 1" 
m . 



(3.11) 
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Let us return to the map induced by the physical flow UV = -6(1,2) 6(2,1) = IR- We 
can write the UV finitized character as 



-—+— 

q 24 ■ 10 



ii,n 

E 



±m uv Cm uv . 



cr,mi,m¥ v 



q \mF v Cm uv 



+ q 2 



\{mi—m\ TV ') 



■|(mi- ) 



9 -I(mi-m^) 







mi — 1 


2 




2 


mi — 1 







(3.12) 



— 1 

2 



+ 



2 



where the three terms correspond to the three mechanisms A, B, C respectively. This is 
confirmed by simple counting in the limit q — > 1 when the Gaussian polynomials reduce to 
binomials. Here mi 



rn 



uv 



m\ does not change under the flow. Assuming (J3T 



and mechanisms A, B, C hold for all N we can uniquely map, q 



q E ' R , energy level by 



energy level. In effect, this means we replace m^ y with ml^ using ( |3~9l) and the energies at 
the base of each Gaussian tower change according to 



A 
B 

C 



q4"- 2 V"" "'2 

q l m uv Cm uv + i {mi+m uv ) 
\m vv Cm uv +|(mi-mf) 



q 4 
q 4 



\rn, IR C™ IR -\ n mi 

q 



2q- 



(3.13) 
(3.14) 
(3.15) 



Putting all this together we obtain the desired mapping between finitized characters 

E 



X(2)(?) 



24 1 16 2 



q \m IR Cm IR 



mi ,m2 R 



r 


N+m! ! R -2 




mi — 1 




2 




2 




mi — 1 




m^ R -l 



+ q 



mi 



N+m^ R - 
2 



2 

mi 



(3.16) 



g 24 + 16 2 



24 1 16 2 



E 

mi ,m2 R 

E 

mi ,77X2^ 




X(2,l)(?) 



7 



4 Boundary TBA Equations B(\2) ^(2,1) 



The derivation of TBA equations is a straightforward extension of OPW to include the 
^-dependent boundary term (|2.4f) . The TBA equations are 



mi 



(1) _ 

ei (a;) = -log^(a;,0-^logtanh(L^)-ir*log(l- e - e2 ^; 

3=1 



tri'z 



t 2 (x) = 4e x — ^ logtanh( 



(2) 



-K* lofffl - e" ei(x) : 



(4.1) 
(4.2) 



k=i 



where * denotes convolution, the kernel is K(x) = 1/ (27r cosh x) and the auxiliary variables 
y% denote the scaled location of the 1-strings in strip i. Similarly, the scaling energies E(j£) 



arc 



E(0 - — = - lim RE(R) = - 



E 



2e~ 



„(!) 



/ — e x log(l 



(4.3) 



where the limit R — > means that all these computations are performed at the critical 
temperature of the TIM. For the lowest level E(£) interpolates between the conformal weights 
A = 1/10 and 7/16 at the boundary fixed points. 

As in OPW the auxiliary variables are determined by a set of auxiliary equations 
determining the locations of the 1-strings 



e 2 (. 7 (1) 



(2) 
k 



■ 7r > 

2' 
2' 



(1) • 
n) in, 

(2) . 



1-strings in strip 1 
1-strings in strip 2 



where n.jj. are given in the IR by the quantization conditions 



(i) 



2 (mi -j + I) 



n 



(2) 



2(m 2 - k + 4 2) ) + 1 



m 2 , 3 ■ 
mi, k 



1,2,... ,m 1 

,m 2 



1,2, 



(4.4) 
(4.5) 



(4.6) 
(4.7) 



The integers nj^ can change during the flow due to winding of phases. The appropriate form 
in the UV is given in OPW. Observe that the boundary field £ appears explicitly only in the 
function g\. 

The 2-strings do not appear in the TBA equations. However, they are involved in 
mechanisms B and C. To follow the movement of 2-strings along the flow we note that their 
locations zf 1 satisfy the same equations ( |4.4|) , ( |4.5|) with the substitution — > zf 1 and the 
appropriate choice of the quantum numbers. We note that in general, and in particular for 
the TIM, it remains an open problem in this approach to obtain the ^-function flow [1], Q 
that interpolates between the boundary entropies ( p..2|) at the fixed points. 
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Figure 3: The flow of scaling energies E(£) between the -6(1,2) an d -6(2,1) boundary fixed points 
as obtained by the numerical solution of the TBA equations. Only scaling energies with A < 7 are 
shown as in Table 1. The indicated degeneracies in the UV (left) and IR (right) are consistent with 
the degeneracies of the Virasoro characters Xi,2(q) = Xi/io(q) an d X2,i{q) = X7/iq{q) respectively. 
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5 Numerical Analysis of the RG Flow Bn 2) >— >• ^(2,1) 



The TBA equations of the previous section can be solved numerically by an iterative proce- 
dure but there are some subtleties. The process starts with an inital guess, for the pseudoen- 
ergies ej(x) and 1-string locations, close to the UV or IR fixed points. The flow is followed 
by progressively incrementing or decrementing the boundary field £. At each value of £, the 
TBA equations are used to update the pseudoenergies e;(x) and then these are used in the 
auxiliary equations to update the locations of the 1-strings, and so on, until a stable solution 
is reached. In our preliminary numerical study we have considered all the states with A < 7 
as shown in Figure [| The corresponding mappings of UV and IR quantum numbers and the 
mechanisms responsible are listed in Table [l|. These flows, and in particular the behaviour 
of the 1- and 2-strings, confirm the mechanisms A, B, C. 

The most important new feature of the numerics is the existence, due to the term fl2.4j) , 
of a pole at x = — £ which marches through the center of the analyticity strip 1 as £ is varied. 
This leads to some inherent numerical instabilities as the pole successively encounters the 
1-string locations in strip 1. A second numerical problem is related to the determination 
of the location of the 1-strings in strip 2. This problem was previously encountered in || 

(2) 

and arises because y k cannot be obtained by direct iteration of the auxiliary equation. It 
is partially solved here using two different schemes 

1. Inverting certain phases containing yp — 

(2) 

2. Inverting the boundary term that contains y k — £. 

The first method, used with success in ||, can be sensitive to errors in the location of y^ . 

The second method ties yjj® to £ so it is generally less sensitive to errors in yj. Usually 
both schemes converge to the same results but in some regions we find one scheme converges 
where the other fails. 
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